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Abstract 

Gravitational physics of VLBI experiment conducted on September 8, 2002 and 
dedicated to measure the speed of gravity (a fundamental constant in the Einstein 
equations) is treated in the first post-Newtonian approximation. Explicit speed- 
of-gravity parameterization is introduced to the Einstein equations to single out 
the retardation effect caused by the finite speed of gravity in the relativistic time 
delay of light, passing through the variable gravitational field of the solar system. 
The speed-of-gravity 1.5 post-Newtonian correction to the Shapiro time delay is 
derived and compared with our previous result obtained by making use of the post- 
Minkowskian approximation. We confirm that the 1.5 post-Newtonian correction to 
the Shapiro delay depends on the speed of gravity Cg that is a directly measurable 
parameter in the VLBI experiment. 
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1 Introduction 



The gravitational VLBI experiment for measuring the relativistic effect of 
propagation of gravity field of orbiting Jupiter was conducted on September 
8, 2002 by the National Radio Astronomical Observatory (USA) and the Max 
Plank Institute for Radio Astronomy (Germany). The idea of the experiment 
was proposed by Kopeikin [1] who had noted that an arbitrary-moving gravi- 
tating body deflects light not instantaneously but with retardation caused by 
the finite speed of gravity propagating from the body to the light ray. This 
is because positions of the gravitating bodies are connected to the light-ray 
particle by the gravity null cone which is defined by the equation for the re- 
tarded time in the gravitational Lienard-Wiechert potentials being solutions 
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of the Einstein equations. The goal of the experiment was to confirm this 
prediction by making use of the close celestial alignment of Jupiter and the 
quasar J0842+1835. The goal of this paper is to make use of the parameterized 
post-Newtonian technique to reveal how the gravity propagation affects the 
light-ray trajectory. This approach allows us to confirm our previous theoret- 
ical results [1] by looking at the problem from different position. 

The problem is formulated as follows. Light ray is emitted by a quasar (QSO 
J0842+1835) at the time to and moves to the network of VLBI stations located 
on the Earth. As light moves, it passes through the variable gravitational field 
of the solar system, created by the orbital motion of Sun and other massive 
planets, and is received by the first and second VLBI stations at the times 
ti and t2 respectively. Gravitational field of the solar system causes delay 
in the time of propagation of radio signals - the effect discovered by [2] in 
static gravitational field approximation. We noted [1] that the present-day 
accuracy of phase-reference VLBI measurements is sufficient to detect a 1.5 
post-Newtonian ^ correction to the Shapiro time delay incorporated implicitly 
into position of the light-ray deflecting body (see Fig.B.l) which must be 
taken at the retarded instant of time related to the time of observation by 
the equation of the gravity null-cone. We solved the Einstein-Maxwell system 
of differential equations and calculated this correction for Jupiter by making 
use of the post-Minkowskian approximation [1] with the retarded Lienard- 
Wiechert solutions of the Einstein equations. The retardation is caused by the 
finite speed of propagation of gravity Cg = c which, as we have revealed in [1], 
enters the 1.5 post-Newtonian correction to the Shapiro time delay. Hence, we 
concluded that the VLBI experiment on September 8, 2002 is sensitive to the 
effect of the propagation of gravity and can be used to measure its speed. 

Our primary paper [1] has given no particular details on the parameteriza- 
tion of the gravity propagation effect the amplitude of which was character- 
ized by the parameter 6. For this reason the physical meaning of the gravity 
propagation effect was not understood adequately by some researchers [3,4]. 
By scrutiny inspection we found that [3] made an erroneous use of insuffi- 
cient 1.0 post-Newtonian (static) approximation to discuss the high-order 1.5 
post-Newtonian (dynamic) correction to the Shapiro time delay while [4] gave 
analysis of the problem of propagation of light ray in the time-dependent grav- 
itational field of the solar system on the basis of the PPN formalism which 
was not sufficiently elaborated to tackle properly the 1.5 post-Newtonian ap- 
proximation to general relativity (see appendix B). 

In this paper we extend our theoretical analysis of the VLBI experiment by 
proving that the parameter 6 = c/cg — 1, where c is the speed of light and Cg 



^ The first 1.0 post-Newtonian approximation is of order c ^, and 1.5 post- 
Newtonian correction to the Newtonian gravity is of order c~^. 
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is the speed of gravity. Parameter 5 is a measurable quantity and its determi- 
nation will provided us with the numerical estimate of the magnitude of the 
propagation of gravity effect in terms of the speed of gravity Cg. The spectrum 
of plausible values of Cg ranges from = c in general relativity to = oo 
as advocated by [5]. The VLBI experiment gives Cg/c = 1.06 ± 0.21 [6]. The 
values Cg < c however are unlikely and ruled out, for example, by observations 
of cosmic rays [7]. 

The paper is organized as follows. In section 2 we discuss the linearized Ein- 
stein equations and the speed-of-gravity parameterization that is used in sec- 
tion 3 for derivation of the post-Newtonian metric tensor. Section 4 deals 
with the post-Newtonian time delay formula for light propagation. The effect 
of retardation of gravity in the gravitational time delay is explained in section 
5. The differential VLBI time delay is derived in section 6 where we prove 
that our previous post-Minkowskian formulation of the VLBI delay is identi- 
cal with the post-Newtonian presentation. Section 7 contains discussion of the 
basic results of the present paper. 

Roman indices run from 1 to 3 while Greek indices run from to 3. Repeated 
indices assume the Einstein summation rule. Indices are raised and lowered 
with the Minkowski metric rjap = diag(— 1, -|-1, -|-1, +1). The Kroneker symbol 
6^ = diag(l, 1, 1, 1). We denote coordinates = (ct, x*) and partial deriva- 
tives da = {c~^d/dt,d/dx^). Partial derivative associated with Cg is denoted 
as Sq = {Cg^d/dr, d/dx^), where r relates to time t via Eq. (3). 

We use the boldface italic letters to denote spatial vectors, e.g, a = a* = 
(a^, a^, a^). The dot between two spatial vectors is the Euclidean scalar prod- 
uct, a-b = a^b^ + a^b"^ + a^b^. The cross between two spatial vectors denotes 
the Euclidean vector product, a x b = Sij^a^b^, where Sijk is the Levi-Civita 
symbol such that £123 = -|-1. 



2 Linearized Einstein Equations and Their Parameterization 



We make use of the following decomposition of the metric tensor g^^ 



where g = det 



and 7^^^ is the gravitational field perturbation. When 



Cg = c the hnearized Einstein equations in harmonic gauge, d^"^"^ = 0, read 



1 



levrG 



x] 



(2) 
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where = d^di is the Laplace operator, and T°'^{t,x) is the stress-energy 
tensor of matter, and G is the universal gravitational constant. 

Special relativity postulates that in every inertial frame, there is a finite uni- 
versal limiting speed c. Experimentally, the limiting speed c is equal to the 
speed of light q in vacuum with amazing precision [8] . General relativity is a 
theory of gravitational field propagating with some speed Cg that was equated 
by Einstein to the fundamental speed c which appear in the wave operator in 
Eq. (2) while the right side of these equations contains a constant IGvrG/c^ de- 
scribing the coupling of gravity with matter. In what follows we intend to build 
up the parameterized version of the linearized Einstein equations (2) with the 
speed of gravity Cg in such a way that allows to get both the Newtonian and 
special relativity limits simultaneously, i.e., — ^ oo while c = const. We em- 
phasize, however, that our approach is just the first step in this direction and 
further development is required. 

Because in general relativity the speed of gravity Cg is denoted by the same 
letter c it is difficult (but possible) to discern the effects associated with Cg from 
those associated with the fundamental speed c. Thus, a theoretical framework 
is needed to unravel Cg from c in the experimental gravitational physics. We 
shall develop in this paper a parameterized approximation to general relativity 
which allows to segregate Cg from c in the experiments conducted in weak 
gravitational fields. Development more advanced (non-linear) formalism goes 
beyond the scope of this paper. Such formalism will be required, for example, 
for setting experimental limits on the speed of gravity from binary pulsar 
timing [9]. 

In what follows, we shall keep notation for the speed of gravity and light as Cg 
and c respectively, and consider the speed of gravity Cg cLS cl parameter which 
can range from c to oo. Speed of gravity characterizes temporal variation of 
gravitational field and the pace of its propagation. This implies that partial 
time derivatives of the metric tensor perturbation 7"^ in the left side of the 
Einstein equations (2) must be parametrically associated with the speed Cg. 
This will lead to the wave operator depending explicitly on the speed of gravity 
Ca considered parameter. 

However, if one holds the time coordinate t fixed for different values of the 
parameter 7^ c it would change the numerical value of the time coordinate 
of one and the same event and mix up the order of space-time points. For 
this reason, we will keep numerical value of the time coordinate = ct to be 
the same while changing c — Cg. This assumes that the smooth transition from 
c to Cg in the left part of Eqs. (2) is to be accompanied with the introduction 
of a new time variable, r, which must scale simultaneously with the parameter 
Cg in such a way that equality CgT = ct is preserved. 
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In other words, the variable r must relate to time t as follows 



r = et , 



(3) 



where e = c/cg = 1 + 6 is a. dimensionless parameter that provides the speed- 
of-gravity parameterization of Eqs. (2) and is used in their post-Newtonian 
expansion to keep track of the orders in the post-Newtonian expansion. As one 
will see later, the time r parameterizes world lines of massive bodies comprising 
a gravitating system. Rescaling of the time variable means that the system 
with e = 1 is doing at time interval AT what the system for any other e 7^ 1 is 
doing at time AT/e, i.e. the faster gravity propagates (e 0) the shorter the 
time scale r. In other words, we are approaching to the Newtonian limit along 
the sequence of solutions of Eqs. (2) parameterized by the sliding parameter e 
but fixed value of = ct. The Newtonian limit is given by e = 0, and in the 
limit e = 1 we restore the linearized Einstein equations (2). 

The left side of Eqs. (2) is parameterized by the speed of gravity after making 
use of the replacements: c — c^, t — > r, and 'j°'f^(t,x) — >• 7"'^(r, a?). This is, 
however, not enough because, if one wants to keep the form of the Eq. (2) be 
formally the same as the original Einstein equations in the harmonic coordi- 
nates, the gauge conditions must be re-written as SQ7°^(r, as) = 0. Since the 
gauge conditions are directly related to the equations of motion of gravitating 
matter [10], their replacement implies that the stress-energy tensor must be 
replaced as well: T"^(t, a:;) 6"^(r, x), where 



This law of re-scaling of the stress-energy tensor is chosen by the requirement 
that the Newtonian limit is obtained when e ^ while the speed of light c 
is constant. As for the formal structure of T°'^{t, x) we assume that it is the 
same as in special relativity. Equations of motion for 6"^(r, x) follow from 
the gauge conditions and are given by S^O"^ = 0. 

Finally, the parametric form of the linearized equations (2) in the parameter- 
ized harmonic gauge, Sa7°^(r, x) = 0, is 



We emphasize once again that it would be illogical to obtain Eq. (5) from the 
Einstein equations (2) with a simple replacement c — > without changing the 
time coordinate t t. An attempt to leave the time coordinate unchanged 
would violate the order of events for different values of the sliding parameter 
Cg and lead to different treatment of the problem of light propagation (see 
appendix B). 
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Smooth transformation of the hnearized equations (2) of general relativity to 
the sequence of parametric space-times allows us to keep two speeds, c and Cg, 
to be formally different for each value of e 7^ 1. The space-time parameteriza- 
tion is introduced via Eqs. (3)-(5) which can be understood as mapping of the 
physical metric gapit.x) (e = 1) into the sequence of solutions QapiT = €t,x) 
of the linearized equations (5), where e is a parameter along the sequence. This 
sequence of solutions is generated for each value of e by the parametric stress- 
energy tensor of matter 0"^(r, x) which is locally conserved. The mapping 
rule for connecting of identical events, having the same spatial coordinates on 
different slices of the parametric space-time manifold, is CgT = ct. 



3 The Post-Newtonian Metric Tensor 



Let us work in the reference frame of the solar system with the origin located 
at its barycenter. Let T""^ be the stress-energy tensor of massive point-like 
particles [10] 



N 

T^^(r, x) = Y: MJ-\t) <(r) <(r) 6^'\x - x„(r)) , (6) 

a=l 

where the index a = 1,2, N enumerates gravitating bodies of the solar sys- 
tem. Ma is the (constant) rest mass of the ath body, Xa{T) are time-dependent 
spatial coordinates of the ath body, Vair) = dXa{T)/dT is the orbital ve- 
locity of the ath body, = 2a{c, Va) is the four- velocity of the ath body, 

2a = (1 — ^ , and 5^^\x) is a 3-dimentional Dirac's delta-function. 

Solution of Eq. (5) can be found in terms of the retarded Lienard- Wiechert 
tensor potentials [11], [12] 



^ ^' ' C r„W-c-i„„W.r„(s)' 

which depend on the retarded time s = s{t, x) determined for each body as a 
solution of the gravity null-cone equation [1] 

(8) 

connecting position of ath body to the field point (r, a:;), and ra{s) = x—Xa{s), 
^ais) = \ra{s)\. Notice that the retarded time s is different for each body as 
it depends on the body's position Xa- 
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The post-Newtonian form of the metric tensor, used in the subsequent calcu- 
lation of a light-ray propagation, is obtained after making use of the Taylor 
expansion of all functions in the left side of Eq. (7) about the present instant 
of time r. It yields 

N ]\/f 

7"(r,-)=4E|i^ + 0(c-) , (9) 

a=l ^ 

7°^(r,x)=46E^<(r)+0(c-^) , (10) 

a=l ^ -"-a 

7'^ (r,x)=4e^5:-— ^<(rK(r) + 0(c-^) , (11) 

a=l ^ ^a. 

where Ra = \Ra\, Ra = Ra = — a;Jj(r). Let us note that in the Newtonian 
limit Cg — s> oo (e — > 0) so that the post-Newtonian metric tensor (9)-(ll) is 
reduced to the simple case with 7°° as the only significant gravitational field 
perturbation. This confirms the standard point of view that gravity propa- 
gates instanteneously in the Newtonian theory and is described by a single 
'scalar' function. It is worth noticing that 7'^'^ directly relates to the metric 
perturbations h^,^ = g^^ - r/^^ = -f^^{T,x) - (l/2)?7^^7"^(r, a;). In the limit 
e = one has /iqo = l/27oo and hij = l/26ij'-fQQ. Thus, integration of the light 
ray propagation equation in this limit will yield exactly the Shapiro time delay 
for the gravitational field propagating instanteneously with infinite speed. 



Calculation of the post-Newtonian limit (9)-(ll) indicates that our parame- 
terized model of the gravitational equations (5) can be interpreted as having 
the following values of the PPN parameters: /3 = 7 = 1, cti = 85, 0:2 = 5(5+2), 
and all other PPN parameters are zero. We emphasize however that this in- 
terpretation makes sense only if one considers the time variable r as physical. 
However, in our formalism time t is physical and coincides with r for e = 1. 
It should be also understood that the velocity Vai^r) = dXa{T)/dT of the grav- 
itating bodies is numerically equal to the physical velocity Vaif) = dXa(t)/dt 
in the approximation of uniform and rectilinear motion which we shall use 
later for analysis of the time delay ^ . This remark is important to avoid mis- 
interpretation of the velocity Vair) in Eq. (23). 

The parameterized post-Newtonian (PPN) formalism [13] predicts that certain 
preferred frame (Lorentz-invariance violating) effects, which are proportional 
to the ratio c/cg, would arise in case if Cg 7^ c. These effects are characterized 
by the product of either of two parameters ai, 02 with the speed w of the 

^ Indeed, if one considers Xa{T) = a + 6r with constant a and b, then Xa{t) = a + ht 
and for any e 7^ 1 one has Va{T) = Va{t) = b. In other words, it is meaningless to 
define velocity as dxa/dt or dxa/dr without explicit indication of the argument of 
the function Xa- 
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solar system with respect to the preferred frame. Under assumption that the 
isotropy of the cosmic microwave background radiation (CMBR) defines the 
preferred frame, a hmit has been estabhshed on ai < 10~^ [13]. But, we do not 
know if the CMBR does, in fact, define the preferred coordinate system. For 
example, future observations of relic gravitational wave background (GWB) 
could lead to the frame moving somehow differently with respect to the CMBR 
frame, since it is assumed that the GWB was formed in the very early uni- 
verse, long before the CMBR decoupled from matter. Furthermore, modern 
multiconnected cosmologies challenge the cosmological Copernican Principle 
and indicates various possibilities for the preferred frame be not coinciding 
with the CMBR frame [14]. We conclude that the existing limit on ai is un- 
certain until we can confidently determine the preferred reference frame and 
the speed w of the solar system with respect to it from laboratory experiments 
measuring violation of the local Lorentz-invariance. 



4 The Post-Newtonian Time Delay 

In general relativity light propagates along null geodesies of the metric tensor 
gapiT, x). In our parameterized approach to general relativity we want to keep 
the speed of light equal to c for any value of the parameter e. This assumes that 
the equations of light propagation must be ordinary null geodesies expressed 
in terms of the Christoffel symbols approximated as 

r^. = Iv''^ {duhp, + d,hp, - dph,,) , (12) 

where /i^^ = ■j^u{,t,x) — ( 1/2)77^,^7"^ (r, a;), and the partial derivative da = 
{c~^d/dt,d/dx^) does not depend on e. The speed of gravity Cg appears in 
equations of hght geodesies and, consequently in the gravitational time delay, 
explicitly through the numerical factors in the components of the metric tensor 
perturbations, and implicitly, through the dependence of coordinates Xa{T) of 
gravitating bodies on time t = et = {c/cg)t. 

Because the speed of light is held constant we must use the coordinate time 
t to parameterize light-ray trajectory while the time r is used for parameter- 
ization of the world-lines of gravitating bodies deflecting the light ray. The 
undisturbed propagation of light in the absence of gravity is a straight line 

x%{t) =xi + ck\t - to) , (13) 

where to is time of emission, Xq is position of the source of light at time to, 
is the unit vector in the direction of propagation of light from the source of 
the light to observer. 
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The post-Newtonian equation of light geodesic parameterized by coordinate 
time t reads 







" ■" ■" V'" ^ ^A'^y ' ^^^^ 

where k^ = is the four-vector of the unperturbed hght-ray path, the 

Christoffel symbols are given by Eq. (12) and should be taken on the unper- 
turbed light-ray trajectory (13). 

Double integration of Eq. (14) yields time of propagation of light from the 
point Xq to x\ = x*(ti) 



1, , A/ N 

ti-to = -\xi - Xo\ + A(ti,to) 



(15) 



where 



A(ti,to 



k^^k^k' 



dt 



to 



dc 



x=XNiO 



(16) 



We emphasize that in order to integrate the right side of Eq. (16) one, first, 
has to take the derivative and, then, to substitute the unperturbed light ray 
trajectory for x. 



One can prove that 



ck'' 



dah^^{et,x) 



d_ 

dt 



{et,XN{t)) 



(17) 



where d/dt denotes the total time derivative along the unperturbed light ray. 
By making use of Eq. (17) one can perform one integration in Eq. (16) and 
recasts it to a simpler form 



A(ti,to) 



-k^'k" 
2 



to 



/i^^ {et,XN{t)) dt . 



(18) 



This equation must be solved for e 7^ 1. It will show at what terms the speed 
of gravity appears explicitly. Then, we take a general relativistic limit e ^ 1. 
This is a way to distinguish the effects associated with the speed of gravity 
(dynamics of gravitational field) from those related to the process of light 
propagation. 

Relativistic time delay is found by performing integration in Eq. (18) with 
positions of the gravitating bodies in the metric tensor (9)-(ll) treated as 
functions of time a?a(r) = Xa{et). The integral (18), however, is too compli- 
cated if the true orbital motion of gravitating bodies is used and its exact 
analytic calculation can not be done explicitly. For this reason, we approxi- 
mate the integral by making use of a linear Taylor expansion 
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x'Ar)=x\{Ta) + vi,{Ta){T - + O (r - r^f 



(19) 



<W=<(^a) + 0(r-r,) , 



(20) 



where t = et and is an arbitrary instant of time that can not be unambi- 
giously determined in the post-Newtonian calculations. In order to fix one 
needs to resort to the post-Minkowskian result of integration of the light-ray 
propagation equation obtained in [1]. This is a very subtle point closely related 
to the interpretation of the relativistic VLBI experiment. It is intuitievly clear 
that Ta should be close to the time of observation but its exact value can be 
found only in the framework of the exact post-Minkowskian theory of integra- 
tion of equations of light, propagating in variable gravitational fields [1]. We 
shall show how to fix the value of in the next section. 

We also emphasize that Xa{Ta) and Va{Ta) are coordinate and velocity of the 
ath gravitating body taken at the time Ta on its real trajectory. Hence, if time 
r in the expansions (19), (20) is significantly different from Tq, the position 
Xair) of the body approximated on the basis of these expansions can be very 
far away from its real position at the time r. This indicates to the limitation 
of the post-Newtonian integration technique which can be cured only within 
the post-Minkowskian approach [11], [12]. 

Substituting Eqs. (13), (19) and (20) for x, Xa^r) and Vai^r) in the metric 
perturbations (9)-(ll) and integrating, we obtain the relativistic time delay 



where Pa = \Pa\, Pa = k-Va/Cg ^ , ilia = Xi - Xai^Ti), Ria = \Rla\, n = eti, 

xi = x(ti), Xa{Ti) and Va must be understood as given by Eqs. (19), (20). 

The constant term Ci {k, xq) is the value of the integral in Eq. (18) taken 
at time to when light was emitted by the quasar. It has no real physical 
meaning in the post-Newtonian approach and could violate consistency of the 
calculation for close celestial objects being members of our galaxy. However, 
the quasar in the relativistic VLBI experiment on September 8, 2002 is lying 
extremaly far away, so that the constant Ci is unimportant in the following 
discussion because it can be ruled out in calculation of the differential VLBI 
time delay that is a real observable quantity. Physically meaningful value of 
the constant of integration Ci can be also obtained only by making use of the 
post-Minkowskian approach [1], [11], [12]. 

^ The speed of gravity Cg appears in the expression for the vector because one 
has x{t) — Xa{T) ~ ckt — eVat = cp^t. Prom the physical point of view it means 
that motion of hght and the bodies in the gravitating system goes on different time 
scales if Cg / c. 



A(ti,to) = -2E-5r^ 

a=l ^ Fa \ J 



In {paR 



Pa-Rla) + Ci (fc, Xo) 



(21) 
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5 The Retardation of Gravity 



By making use of the post-Minkowskian theory of light propagation we have 
predicted [1] that the VLBI experiment on September 8, 2002 is able to mea- 
sure the speed of gravity by observing the effect of the retardation of gravity 
propagating from the sun-orbiting Jupiter to light ray. In this section we con- 
firm this interpretation on the basis of the post-Newtonian approach. To this 
end, let us recast Eq. (21) to another form by making use of the Taylor ex- 
pansion of vector with respect to the small parameter Va/cg. This yields 



TV 



A(ti, to) = -2 E —J^ 1 - -k-va In 



a=l \ '^9 



Pa {Rla - ria-Rla) 



(22) 



where C (fc, a;o) is (another) constant of integration and the unit vector ^ 

na = k — -k X {va X k) (23) 

'^9 

explicitly describes a change in the direction of propagation of light that can 
be interpreted as dragging of the light-ray caused by the flux of the gravi- 
tational fleld generated by the translational motion of ath body. This phe- 
nomena can be equivalently interpreted as the effect of aberration of gravity, 
i.e. displacement of the observed body from its present to retarded position, 
which is unobservable in the equations of motion of slowly-moving particles 
[15] but can be detected in the light-ray deflection relativistic experiments in 
time-dependent gravitational flelds. 

This can be demonstarted if one chooses the time Ta in Eqs. (19) and (20) 
to be coinciding with the retarded time si corresponding to the instant of 
observation ti and reflecting the fact that gravity takes flnite time to propagate 
from the moving gravitating body to observer. This time is calculated from 
Eq. (8) as follows 

.i = n-^, (24) 

Cg 

where n = eti, ria(si) = Xi - Xa{si), and ria(si) = |ria(si)|. The geometric 
distance Ria between the observer and the ath body reads 



Rla = riaiSl) - ^^i^ Vais,) + 0(0^) , (25) 

Ria = rUsi) - -v,{s^)-rUsi) + 0{c;^) . (26) 



^ We again emphasize that the velocity Va = i'a('^a) is the physical velocity of the 
ath body. 



11 



Substituting Eqs. (25) and (26) to Eq. (22) and reducing similar terms one 
obtains 



A(ti,to) 



a=l ^ 



1 k-Va ) In 



ria{si)-k-ria{si) 



-Cik,xo) . (27) 



This formula gives the post-Newtonian time delay in terms of the retarded 



position of the solar system bodies and coincides in the limit Cg 



with 



the analogous formula derived by [1] on the basis of the post-Minkowskian 
approach ^ . It is evident from Eq. (24) that the retardation in Eq. (27) is due 
to the finite speed of propagation of gravity Cg because Eq. (24) is just the 
equation of the gravity null cone describing propagation of gravity between 
the moving gravitating body and the point of observation. Relativistic VLBI 
experiment on September 8, 2002 measured the retarded position of Jupiter in 
the sky through the gravitational deflection of light of the quasar which was 
in the plane of the sky at the distance of ~ 14 jovian radii from Jupiter (see 
Fig. B.2). 



6 Discussion 



We have constructed a speed-of-gravity parameterization of the linearized 
equations of general relativity and calculated a speed-of-gravity correction 
to the Shapiro time delay given by Eq. (27). Its post-Newtonian expansion 
is identical with Eq. (12) from our previous work [1] where we used the 
post-Minkowskian approach and introduced the fitting parameter 5 on a phe- 
nomenological basis. The present paper proves that on the sequence of the 
parametric space-times manifolds, 6 = c/cg — 1 and distinguishes two funda- 
mental speeds in the Einstein and Maxwell equations. Hence, the relativistic 
VLBI experiment conducted on September 8, 2002 allows to determine the 
speed of gravity and to measure its numerical value. 

Einstein theory of general relativity predicts that perturbation of the gravita- 
tional field produced by orbital motion of Jupiter propagates with the speed 
of gravity Cg = c towards the position of the light-ray particle traveling from 
quasar to Earth. This prediction can be tested by VLBI technique without 
direct detection of gravitational waves. The VLBI experiment gives the magni- 
tude of the retardation-of-gravity effect in terms of its speed Cg/c = 1.06±0.21 
[6]. We conclude that the speed of gravity is the same as the speed of light 
within our experimental error. 



Let us emphasize that if one did not choose ti = si such a coincidence could not 
be achieved. 
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A Comments on the paper by Asada [3] 

Asada [3] calculated the relativistic time delay which formally coincides with 
our formula (27). Asada's calculations are not complete as he postulated in his 
calculations that position of Jupiter must be fixed at the retarded instant of 
time si = ti — (l/c)|a;i — xj{si)\. Such postulate is a plausible assumption 
but not a mathematical proof. For this reason Asada's interpretation of the 
relativistic VLBI experiment on September 8, 2002 is ungrounded. 

In our previous paper [1] we have solved the system of the Einstein and light 
geodesic equations without making any restrictive assumptions about the mo- 
tion of Jupiter and proved that the position of Jupiter must be taken at the 
retarded instant of time Si which is an argument of the metric tensor and, 
hence, describes the effect of propagation of gravity with the gravity speed 

Cg = c. 



B Comments on the paper by Will [4] 

Will [4] calculated the relativistic time delay by making use of the matched 
asymptotic technique of the solutions of the light geodesic equation solved in 
the near- and far-zone fields of the gravitational potentials U and Vi obtained 
from the field equations 




(B.l) 




(B.2) 
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where p(t, x) is the density of matter and the time t is the same as in the hght 
geodesic equations. 

Eqs. (B.l), (B.2) are not Einstein equations and are written in a particular 
gauge which was not clearly described. Furthermore, we notice that the field 
equations (B.l), (B.2) are, in fact, disconnected from the light geodesic equa- 
tions. Indeed, assuming that c and using a single time scale t imply that 
vacuum must be considered as having a "refractive" index n ^ cj Cg ^ \ and 
the equations of light propagation should have extra terms proportional to 
n — 1 [16] that were not discussed in [4]. The whole problem of the unam- 
biguous interpretation of the VLBI experiment under discussion is, in fact, 
a matter of construction of a self-consistent bi-metric theory of gravity and 
electromagnetism. The PPN formalism [13] can not play a role of such a the- 
ory as it consists of a set of plausible but, in any other aspect, disconnected 
assumptions about how gravity and electromagnetic fields behave if general 
theory of relativity is violated. 

There are two other concerns about validity of the calculations given by Will 
[4]. The first concern relates to the meaning of the retarded time = cr — 
faip-, So) /cg in the arguments of the bodies Xa{sa)- Will insists (see appendix of 
[4]) that coefficients of the Taylor expansion Xa{sa) = Xa{S^)+Va{S^)sa + ■■■ must 
be understood as taken at the time a = 0. This contradicts to the definition of 
the Taylor expansion given in any mathematical textbook. By definition of the 
Taylor expansion the coefficients a;a(0) and i'a(O) in the formula given above 
are taken at the retarded time Sq = ^ . This instant of time can coincide 
with that cr = 0, if and only if, the retarded time equation Sa = — ?"a(0, Sa)/cg 
has a solution Sa = 0. From the astronomical point of view it means that 
the light ray must pass through the light-ray defiecting body. This is indeed 
coincides with what is shown in the figure of Will's paper. However, in the 
relativistic VLBI experiment on September 8, 2002 light was passing pretty 
far away from Jupiter at the distance of about 14 jovian radii (see Fig. B.2). 
Hence, for cr = the time Sa 7^ in the VLBI experiment. 

The second concern relates to the definition of the point of the closest approach 
of the light ray to the moving body. After removing the misinterpretation of 
the Taylor expansion in [4] one has to understand that Will has postulated 
that the point of the closest approach is reached by light at the instant of 
the retarded time Sa = 0. Will also postulated that at the point of the clos- 
est approach the light-ray-propagation unit vector k and that of the impact 
parameter da{sa) = $, — Xa{sa) must be orthogonal in the Euclidean sense. 
We emphasize that the point of the closest approach of light ray to the body 



° Misinterpretation of the retarded time arguments led Will to the erroneous use 
of the Taylor expansion in his Eq. (34) that, as a matter of fact, must be used along 
with the gravity null cone equation (18). 
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must be found after solution of the corresponding equation giving the dis- 
tance between the body and the hght particle as a function of time, let say, 
D(cr, Sa) = \cka + da{sa)\- Finding the minimum of this function from the con- 
dition dD{a, So) /da = 0, shows that at the point of the closest approach the 
two vectors k and da{sa) are not orthogonal and their Euclidean dot product 
is proportional to the relativistic terms of order daVa/cg which are important 
in the subsequent calculations of the time delay. All such terms are absent 
in [4] that along with the previous our comments casts strong doubts on the 
validity of Will's calculations and his interpretation of the relativistic VLBI 
experiment conducted by Fomalont and Kopeikin [6]. 
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Light from Quasar 




Fig. B.l. Light propagates with speed c from the source of electromagnetic waves 
(quasar) to radio observatory. A massive body (Jupiter) orbits the Sun with speed 
V. Perturbation of gravity field produced by the orbital motion of Jupiter propagates 
from Jupiter to the hght ray with the speed Cg = c. In particular, the gravity-field 
perturbation arrives to the point of observation of radio signal at time ti . In order to 
calculate the delay in the time of arrival of light from the quasar to observer, position 
of Jupiter must be calculated at the retarded time si = ti — r/cg irrespectively of 
the direction of propagation of light [1]. The present position of Jupiter at the time 
ti is displaced with respect to its retarded position at the distance (f/cg)r cos -0 and 
does not play any role in the calculation of the time delay of light. Sun and Earth 
produce similar effect. 
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Jupiter's Worldline Light from Quasar 



Fig. B.2. Minkowski diagram showing the relative configuration of gravity and light 
null cones related to the point of observation. Two cones do not intersect but are 
tangent to each other along a single (dashed) line. The tangent line of the two null 
cones connects Jupiter and observer. Light from quasar propagates along a different 
world line on the light null cone emanating from the quasar (not shown). The 
gravitational interaction between moving Jupiter and a light particle always occurs 
along the gravity null cone. In other words, when Jupiter moves it continuously 
sends gravitational "messages" about its position in space that propagate along the 
gravity null cone with the finite speed Cg = c towards the light particle. The last 
such "message" reaches the light particle at the point of observation and points 
out to the retarded position of Jupiter. This interplay between gravity and light 
was used in the VLBI experiment on September 8, 2002 to test how fast gravity 
propagates and to measure its speed [1], [6]. Will's statement [4] is that position 
of Jupiter in the VLBI time delay formula (27) must be taken at the time of the 
closest approach of light to Jupiter. 
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